Fast intersystem crossing in transition-metal complexes 
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The mechanism behind fast intersystem crossing in transition-metal complexes is shown to be a 
result of the dephasing of the photoexcited state to the phonon continuum of a different state with 
a significantly different transition metal-ligand distance. The coupling is a result of the spin-orbit 
interaction causing a change in the local moment. Recurrence to the initial state is prevented by the 
damping of the phonon oscillation. The decay time is faster than the oscillation frequency of the 
transition metal-ligand stretch mode, in agreement with experiment. For energies above the region 
where the strongest coupling occurs, a slower "leakage" -type decay is observed. If the photoexcited 
state is lower in energy than the state it couples to, then there is no decay. 

PACS numbers: 78.47.J-, 33.50.-j 



Introduction.— Fast intersystem crossing is an in- 
triguing phenomenon that has puzzled many for several 
decades [l|-l3l • In a wide variety of transition-metal com- 
plexes, laser excitation creates a photoinduced excited 
state that decays on the order of tens to hundreds of 
femtoseconds into a state with often a different spin and 
a significant change in transition metal-ligand distance. 
This is generally followed by a cascade of intersystem 
crossings. Examples are spin crossover phenomena in 
divalent iron. The iron atom is generally surrounded 
by an organic material, such as Fe[(phen)2 (NCS)2] or 
Fe(bpy)3 (bpy=bipyridine) In the ground state, 

Fe^+ is in a low spin state. Illumination by light causes 
a charge transfer to the ligands, followed by a cascade of 
intersystem crossings, turning the singlet configuration 
{S = 0) into a high-spin S — 2 state. The high-spin 
state has a temperature-dependent decay time back into 
the low-spin state that can vary from nanoseconds up 
to days. Comparable crossovers have been observed in 
nickel compounds, where the electronic transition is fol- 
lowed by dissociation of the molecule Q . These phenom- 
ena have been explored since the sixties predominantly 
using optical and Mossbauer spectroscopy [3]. Recent 
ultrafast x-ray spectroscopic measurements |^-|6[ have 
provided more detailed information on the crucial first 
step that sets up the cascading process showing that the 
electronic transition occurs on the order of tens of fs [1] 
whereas the lattice relaxes on a timescale related to the 
transition metal-ligand stretch mode (100-200 fs) Q. 

Experimental work has focused on identifying the 
states and relevant times scales of the intersystem cross- 
ing. Theory has followed several approaches. Study of 
the time-dependence generally relies on phenomenologi- 
cal rate equations Energy level diagrams have been 

obtained using ab initio techniques that use an adiabatic 
approximation [lol [Tlj , which has the disadvantage that 
strong coupling between states only occurs for very par- 
ticular energies and transition metal-ligand distances. 

This Letter focuses on the important first intersystem 
crossing and explains why the system decays so rapidly 



into states that dramatically differ in both spin and lat- 
tice parameter with quantum efficiencies close to 100%. 
Several types of decay depending on the relative energy 
positions of the states are identified. 

Excitation.— Figure [1] gives a schematic diagram of the 
fast intersystem crossing. Initially, the system is in the 
ground state \ipo)- It is then excited into state \ipi) from 
which it will decay into state |'02)- Depending on the 
change in coupling to the ligands, several phonons can 
be created in the excitation process, leading to states 
{ipin). Although details of the initial excitation can play 
a role, its primary role is to create an excited state higher 
in energy compared to states that have a different spin. 

Dephasing.— We now turn our attention to the crucial 
step of the fast intersystem crossing. Transitions between 
different spin states are accompanied by a change in lat- 
tice constant due to the conversion of t2g into Cg electrons 
that repel the ligands more strongly (for example, from 




FIG. 1: (color online) Schematic of fast intersystem crossing. 
The initial photoexcitation from state \^o) into state \il^in) 
occurs with minimal lattice distortion, i.e., the number of 
excited phonons n is small. State IV'in) couples to state \ip2n) 
under the excitation of phonons. The coupling strength (red 
curve) creates an effective phonon density of states. 
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LS to HS t^gCg). THs can be described by the Hamil- 
tonian 

i?o = ^ Eic\ci + hwa^ a + ^ eihwuila^ + a), (1) 



2=1,2 



i=l,2 



where iJ^ = £'2 + A, £'2 for the relative positions of states 
1 and 2, respectively. The third term describes the cou- 
pling between the electrons and the lattice. Classically, 
this term corresponds to a constant force displacing the 
equilibrium position of the ligand atoms. Only the rel- 
ative change in coupling is of importance. Since exper- 
imentally one observes a change in spin of AS" = 1, the 
coupling, 



Hi=V{c\c2+c\ci). 



(2) 



is generally accepted to be due to the spin-orbit coupling. 

We can diagonalize Hq with a displaced-harmonic- 
oscillator transformation i/o = e^Hoe~^ with S = 
J2i ni^/gl{a^ - a) with = eijhw. 



giving 



Hq = ^(-Bj - ei)clci + fihJo)a, 



(3) 



with eigenstates |V'm) for states i = 1,2 and n ex- 
cited phonon modes. The transformed coupling Hi = 
e^Hie~^ with Vnn' = {'^2n\Hi\ipin') is no longer phonon 
conserving. If the change in transition metal-ligand dis- 
tance in the photoexcitation is small, the prepared state 
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FIG. 2: (color online) The probability Pi (red) of finding state 
1 and the relative change in transition metal-ligand distance 
(blue) with x/xi = 1,0 for the equilibrium positions of 1 and 
2, respectively, as a function of the time t after the photoex- 
citation. (a) No phonon damping; Sij indicates the change 
in coupling to the lattice between states i and j. Here the 
change in lattice occurs between states 1 and 2; (b) idem, but 
with phonon damping; (c) with phonon damping, but now 
the strongest change in lattice parameters occurs in the pho- 
toexcitation from to 1 (only Pi is shown). The total change 
is £02 = 0.2 eV and £oi = 0.2 (solid) and 0.15 (dotted) eV. 



is close to IV'io), see Fig. [T] The dominant coupling to 
\tp2n) is then V^^ = e'l^g'^/nl with ^ = ^ - 



The strongest coupling between \ipio) and \1Jj2n) occurs 
at the maximum of the Poisson distribution for which 
n = g at an energy E2 + ghuj ~ E2 + £12. In the calcu- 
lations, we take hu = 30 meV, which is a typical value 
for a transition metal-ligand stretching mode The 
value of £12 = £2 — £1 = 0.2 eV corresponds to a typi- 
cal displacement of a few tenths of an Angstrom [ll| . 
This energy difference also corresponds well to the ob- 
served change in wavelength of the emission features in 
the luminescence spectra [1]. The spin-orbit coupling 
parameter is ^ = 50 meV, which is the atomic value 
calculated in the Hartree-Fock limit (T^]. Let us first 
consider A = 2£i2, when ji/'io) coincides with this maxi- 
mum, see Fig. [T] FigurelUa) shows the time-dependence 
of the probability Pi of finding the system in state 1, with 
= l(V'm|V'(^))Pj and the displacement of ligands 
{x ~ {ip{t)\x\tp{t)), where x/xi=l, correspond to the 
equilibrium positions of states 1 and 2, respectively). In 
the first 20 fs, we observe a rapid decrease in Pi due 
to the dephasing of \tpio) into the phonon states |i/'2n)- 
Although dephasing is more commonly associated with 
coupling to a continuum [l3| , there are sufficient phonon 
states for a fast dephasing. However, the finite spacing 
huj between the phonon states leads to a recurrence as 
evidenced by the increase in Pi around 140-160 fs. This 
recurrence is directly related to the oscillation of the lig- 
ands, where hu — 30 meV corresponds to a period of 138 
fs. When looking at larger times, we see a finite value of 
Pi , which does not correspond to the experimental obser- 
vation of quantum efficiencies close to 100%. In addition, 
there is a continued oscillation of the ligands. 

Damping.— The quantum efficiency can be dramati- 
cally improved by including damping of the phonon os- 
cillations [15| . Here we describe how to include a damp- 



ing to a bath that accounts for dam pin g mainly due to 
intramolecular energy redistribution |l6l . flTj . In the ab- 
sence of the coupling Hi relaxes the oscillations into their 
equilibrium states lipio). Since fast intersystem crossings 
are experimentally known to be almost temperature in- 
dependent, we take T = 0. The wavefunctions \ip{t)) = 
J2inb ^ini^)\'^in) ' ^^'^ extended to include the possible 
states b of the bath. The bath variables are integrated out 
to obtain the time-development of the occupation of state 

\i'^n), P^n{t) = E J (V'L IV' W) P = EJcJ„(i)l'- For the 

matrix elements of the coupling between the local states 



and the bath, (V'^,„_i|i?oB|V'i 
the terms linear in a and 



\/nVb5ij, we take only 
The superscript b in the 



states IV'in) indicates that an additional excitation with 
energy Ef, has been created with respect to state with the 
original bath (I'f/'in))- From the Schrodinger equation, we 
find for the change in coefficients due to presence of the 
bath 



dCin{t) 



dt 



(4) 



(for clarity, we omit the subscript B in the following equa- 
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tions). We omit processes that increase the number of 
phonons n via processes involving and the creation of 
an excitation in the bath (the rotating-wave approxima- 
tion). These processes have an energy hu+Et, and are less 
likely than the phonon conserving excitation of energy 
huj—Eb. Numerical calculations [18] show that the proba- 
bility of these processes are an order of magnitude smaller 
for the damping time scales under consideration in this 
paper. The coefficients on the right-hand side of Eqn. 



dl]), can be found from dci^_^{t)/dt 



.JnVbC^n{t). 



Integration gives for the change in the norm Eqn. ^ 



d\cin{t)\ ^ 
dt ^ 

b 



x-l(t-*'). 



At'), 



where = Ei — Si + Ei, is an eigenstate of Hq and 

the bath Hamiltonian. The integral is greatly simplified 
when using the common approximation that the sum over 
the bath integral can be replaced by an integral over an 
effective density of states. Taking a density of states p 
and coupling constant V in the region where it couples 
strongly, we obtain 



1/2 

b 



*') 



ttpbV^ 



ft2 



S{t-t'). 



(5) 



Introducing the lifetime broadening F — irpBv/h, we 
obtain dcin/dt ~ —nTcin{t) . This term shows the de- 
crease in Cin with the bath unchanged in the time step. 
However, at the same time Ci„ can increase by the decay 
of the state with n + 1 excited phonons. By calculating 
the probability for „ and using Eqn. ([5]) , we obtain 



dt 



= -2nrP,r.{t) + 2{n + l)TP,^n+i{t) (6) 



The obtained rate equation can also be obtained heuris- 
tically by taking an oscillator frequency with damping 
U} — ll hjl. The chan ges in the coefRcients due to damp- 
ing is included in the Schrodinger equation through [18| 



ih 



d\ m) 

dt 



ihdlnP{t) 
~2 di 



\m), (7) 



with p{t) = Ez„^»WIV'»«)(V'»«|. 

Experimentally, no oscillations of the ligands are ob- 
served indicating that the intramolecular vibrational en- 
ergy distribution occurs within half an oscillation period 
of the stretching mode. Figure [Hb) shows that the in- 
clusion of F^^ = 30 fs causes a damping of the oscilla- 
tion within a few periods (for the rest, using the same 
parameters as Fig. IHb))- The damping strongly su- 
presses the recurrence of intensity in state 1. The net 
result is a very fast intersystem crossing with a close to 
100% quantum efficiency in good agreement with exper- 
iment. The decay time, however, is not simply related 
to the oscillation period of 138 fs of the phonon, but sig- 
nificantly smaller. A good estimate of the decay time 
is obtained by viewing it as a dephasing of state |'0io) 



into the phonon continuum of that |'02n)- The coupling 
strengths Vno = e~^g"-/n\, can also be viewed as a con- 
stant coupling F to a density of states, see Fig. [U given 
by p{e) — e~'^g'^ /[hujr{e + 1)], with e ~ E/hu! and where 
F(e -I- 1) is the gamma function that allows calculation of 
the factorial for non-integer values. The lifetime broad- 
ening Td associated with the dephasing is then straight- 
forwardly obtained [14] through F^ = ttV'^p{A — e)/h. 
This gives a decay time of (2Fc;)~^ = 26 fs. The ex- 
ponential decay e~'^^'^* provides a satisfactory agreement 
with the numerical calculation, see the dotted line in Fig. 
[2Ib). Furthermore, this rapid decay does not require the 
involvement of high-frequency phonon modes (150-200 
meV), as was suggested earlier since the electronic and 
lattice relaxation occurs on different time scales. 

It has also been suggested 3 that the major change in 
transition metal-ligand distance occurs in the photocxci- 
tation. Figure [2jc) shows the calculation for that situa- 
tion under the assumption that the total change is con- 
stant (eo2 = 0.2 eV). If there is no difference in bond 
length between states 1 and 2 (ei2 — and £02 — 0.2), 
one observes oscillations in Pi. Note that in this limit, 
phonons are created in the photoexcitation. For £12 = 0, 
there is no phonon continuum but a simple coupling be- 
tween the two states. When introducing a small differ- 
ence in lattice between states 1 and 2 (£12 = 0.05, see 
dashed line in Fig. HJc)), we do not find a fast decay nor 
a close to 100% quantum efficiency. 

We now turn our attention to the dependence of the 
relative positions of states 1 and 2. When using an adi- 
abatic approach, strong coupling between the different 
states only occurs for very particular energies and tran- 
sition metal-ligand distances. This appears in contradic- 
tion with the prevalence of fast intersystem crossings. We 
demonstrate that the fast decay occurs for a relative wide 
energy range. In Fig. [31 we show the decay and change 
in transition metal-ligand distance for several values of 
A. We take that the initially photo-prepared state is 
I'f/'io), i.e. the change in the lattice upon photoexcitation 
is small. As described above, the fastest decay occurs for 
A = 2£i2, see curve e in Fig. |3l However, also in a re- 
gion more than £12 wide around this point (see d and f), 
we see the same dephasing with the lifetime broadening 
following the effective density of phonon states p{E), see 
the lower part of Fig. [3] When going to curve g, there is 
an initial dephasing followed by a slower decrease in Pi. 
However, even at higher energies (h and i), where the ef- 
fective phonon density of states p{E) is so small that de- 
phasing does not occur, the decay occurs on a timescale 
larger than the oscillation period of the phonons. We 
can understand this behavior as a "leakage" into state 
2. Since the excitation is into \ipio), phonon damping is 
irrelevant for state 1. The strongest coupling to state 2 
occurs at the maximum of the effective phonon density 
of states. For significantly large A, we can treat p{E) as 
a single state and write an effective Green's function for 
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FIG. 3: (color online) The probability Pi (alternating dark 
blue and dark red) of finding state I and the relative change 
in transition metal-ligand distance (alternating dotted blue 
and red) xlx\ with the f and values indicated on the left 
and right axes, respectively, in corresponding colors. The de- 
pendence as a function of time t is given for different positions 
of I'f/'io) with respect to the effective phonon density of states 
p[E) of state 2, see bottom part of the figure. 
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(8) 



with z = E— E2 — iO"*". The phonon damping F in- 
directly affects state 1. The effective state in the Green's 
function consist of \1p2n) with n = 812/ huj. Phonon 
damping causes these state to relax into | "020)1 which cou- 
ples very weakly to IV'io)- This causes a decrease in Pi 
with an effective decay time aF^F/[(A — £12)^ + [aVY]. 
The factor a accounts for the fact that there are ap- 
proximately g phonons excited. For an independent 
oscillator, the rate equation can be solved analytically 



and is equivalent to g independently decaying oscilla- 
tors with one phonon excited. The time needed to ob- 
tain a 50% occupation of the n = state is 1/aF with 



= -ln[l-(l/2; 



1/91 



When state iV'io) is on the low-energy side of the ef- 
fective phonon density of states p{E) (a-c in Fig. [S]), we 
see an entirely different behavior. When A = ei2, states 
I'f/'io) and |'02o) are almost degenerate, and one obtains 
a strong coupling between them. For, A > £12, state 

1 slowly decays into state 2. For A < £12, state 1 does 
not decay into state 2. At first, the asymmetry below and 
above the phonon continuum seems surprising since IV'io) 
can still couple to the effective phonon continuum of state 

2 and one would expect leakage. However, now state 2 
can easily couple back into the phonon continuum \il)in) 
of state 1, which is now at a much lower energy (about 
£12 above The states then damp into IV'io)- 

Conclusions.— An explanation is given for the fast in- 
tersystem crossing in transition-metal systems. The large 
change in lattice between the states produces an effec- 
tive phonon continuum leading to a fast dephasing. The 
change in spin is necessary since most parts of the elec- 
tronic Hamiltonian do not produce an intersystem cou- 
pling. Damping of the phonon oscillation by coupling 
with the environment prevents a recurrence of intensity in 
the initially photoexcited state. We demonstrate that an 
intersystem crossing of tens of fs is compatible with a lat- 
tice relaxation of the order of 100-200 fs. Different regions 
of decay can be identified as a function of the relative 
energy positions of the states. The fastest decay is an ef- 
fective dephasing and occurs on a smaller timescale than 
the oscillation period of the phonons that are involved 
in the process. Quantum efficiencies of close to 100% 
are obtained. After obtaining an understanding of the 
crucial initial step leading to a complete spin crossover, 
future research should provide a quantum-mechanical un- 
derstanding of cascading effects and finally close the loop 
with a return to the ground state. Additionally, spectro- 
scopic signatures of the process need to be investigated. 
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